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It is shown that Newtonian mechanics is not appropriate to compute approximate individual trajectories with 
small velocities in chaotic scattering. However, some global properties of the dynamical system, such as the 
dimension of the non-attracting chaotic invariant set, are more robust and the Newtonian approximation provides 
reasonably accurate results for them in slow chaotic scattering. 
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I. INTRODUCTION 

It has been recently pointed out that, when the evolution 
is chaotic, the trajectories predicted by Newtonian mechan- 
ics are not good approximations to the special-relativistic tra- 
jectories, even for very small velocities, fljj— 13|] . The problem 
considered in Refs. UHM is the evolution of periodically delta- 
kicked systems with or without damping. The purpose of this 
letter is twofold: we want to extend the previous conclusion 
to another area of deterministic chaos (chaotic scattering) but 
also to point out that, for slow motion, Newtonian mechanics 
can provide good approximations for global characteristics of 
chaotic systems, such as the structure and dimension of the 
non-attracting chaotic invariant set. 

We will consider a relativistic particle moving in an exter- 
nal potential energy V(r). If particle's mass, position and lin- 
ear momentum are m, r and p, respectively, the Hamiltonian 
(i.e., the total energy) is H = y m 2 c 4 + c 2 p 2 + V(r) and 
Hamilton's canonical equations read as follows: 



Relativistic scattering 



r = — , p = -W(r), 
where the Lorentz factor is given by 
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(1) 
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The Newtonian equations of motion are recovered by putting 
7 = 1, so that one might think they are a good approximation 
for slow motion, where 7 w 1. We will discuss in the next 
section to what extent is right this assumption in the case of 
chaotic scattering. 



H. CHAOTIC SCATTERING 

Let us consider a particle moving in the plane (x, y) under 
the force derived from the potential 10, Ht] 



V(x,y) = x 2 y 2 e-( x2+y2 ). 
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Newtonian approximation 
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FIG. 1: Scattering angle <!>(&) for the relativistic and Newtonian 
equations of motion. 



The particle is sent from (x, y) — (—00, b) with initial ve- 
locity (x,y) = (0,v) and after some time emerges from the 
scattering zone moving freely along a direction forming an 
angle $ with the positive x axis. The initial velocity v is com- 
puted from the conserved energy E = H — mc 2 . We choose 
a system of units in which m = c = 1 and a value of the 
energy E — 0.05Vb — where Vq = e~ 2 is the height of the 
four maxima of the potential located at (x, y) = (±1, ±1) — , 
so that the motion is slow, since the Lorentz factor 7 is always 
smaller than 1.007. 

In Fig. [1] we have depicted 1000 values of the scattering 
function <&(&), in the range [0.2, 0.25] of the impact parame- 
ter b, obtained by numerical integration [6] of the relativistic 
equations of motion (fl} and their Newtonian approximation. 

First of all, one can see that while Newtonian mechan- 
ics gives reasonably approximate results when the scattering 
function is smooth, when the latter is discontinuous the pre- 
dictions of relativistic and Newtonian equations of motions 
are very often completely different. This extends to chaotic 
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scattering the claim in Refs. IH-dtl that Newtonian equation 
of motions are not a good approximation when individual 
chaotic trajectories have to be computed. 

However, although the positions of individual points in the 
regions of chaotic scattering may be very different, the posi- 
tion of the discontinuity points along the b axis seems to be 
very similar in both cases. This statement can be made quan- 
titative if one computes the dimension of the Cantor set of 
discontinuity points by using the method of the uncertainty 
exponent (4j,|5[]. We have computed the scattering function for 
2 22 values of b € [0.2,0.25] and then selected a small value 
e. We say that a value b is certain if for the three impact pa- 
rameters b — e, b and b + e the scattering is upward or if it is 
downward for the three values. We repeat this for many val- 
ues of b and call /(e) to the proportion of uncertain values of 
b. Figure|2]show a log-log plot of several values of /(e) along 
with a good linear fit, so that /(e) ~ e a with a uncertainty 
exponent given by the slope of the straight line: a f=s 0.139. 
In consequence [4, 5], the fractal dimension of the set of dis- 
continuity points is Dq = 1 — a s» 0.861. The corresponding 
plot for the Newtonian approximation looks nearly identical 
with a w 0.141. This shows that the dimension of the chaotic 
invariant set is well approximated when computed with New- 
tonian mechanics, if the velocities are small enough. 
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FIG. 2: Log-log plot of some points of the uncertainty function /(e) 
and a linear fit. 



III. FINAL COMMENTS 

In the previous section we have used Newtonian mechanics 
in a consistent way, so that for each impact parameter b the 
initial velocity v is computed using the Newtonian definition 
for the kinetic energy. Since the latter is different in special 
relativity, slightly different initial velocities were used in the 
Newtonian and the relativistic simulations. If one uses the 
same initial velocities in both cases (which implies slightly 
different energies), the results are very similar to those dis- 
cussed above. For instance, one gets a rs 0.142. 

We have shown that also in chaotic scattering Newtonian 
mechanics fails to be a good approximation to compute indi- 
vidual trajectories. Even for slow motion, the little relativistic 
corrections to the equations of motion, which remain small 
in non-chaotic scattering, are greatly amplified by chaotic 
scattering during the long time in which the trajectory re- 
mains close to the non-attracting invariant set. This extends 
to chaotic scattering the conclusions reached in Refs. JTrH]. 
However, we have also shown that some global properties, 
such as the structure and dimension of the set of discontinuity 
points (and, thus, the dimension of the non-attracting chaotic 
invariant set) are more robust and Newtonian mechanics pro- 
vides reasonably accurate approximations to them if the mo- 
tion is slow enough. 
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